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We reinvestigate a doped antiferromagnetic Mott insulator based on the slave-fermion approach of 
the t-J model, where antiferromagnetic spin fluctuations and doped holes are described by bosonic 
spinous and fermionic holons, respectively. Earlier studies have shown that an effective field theory 
for the doped antiferromagnetic Mott insulator is given by a non-relativistic fermion (holon) U(l) 
gauge theory for charge dynamics and a relativistic boson (spinon) U(l) gauge theory for spin 
dynamics, thus allowing an anomalous metallic phase where bosonic spinous are gapped away from 
an antiferromagnetic state, analogous to the U(l) spin liquid phase in the slave-boson approach 
of the t-J model. We argue that the emergent U(l) gauge structure in this approach is based 
on a simplified picture for antiferromagnetic correlations. Considering that dynamics of doped 
holes frustrates a collinear antiferromagnetic spin configuration, we show that doped holes enhance 
ferromagnetic spin correlations and result in a complex magnetic structure. The presence of such 
a complex spin texture reduces the U(l) gauge structure down to Z2, thus a different effective field 
theory results for the spin-gapped non-Fermi liquid metal at low temperatures because there are no 
gapless U(l) gauge fluctuations. 

PACS numbers: 71.10.Hf, 71.10.Fd, 71.27.+a, 75.10.-b 



I. INTRODUCTION 

Study on strongly correlated electrons opened a new 
window so called gauge theory in modern condensed mat- 
ter physics. When electrons are weakly correlated, that 
is, kinetic energy is more dominant than potential en- 
ergy, a local order parameter approach is available based 
on decoupling of interaction channels. This is the heart of 
" classical" condensed matter physics composed of Fermi 
liquid theory and Landau- Ginzburg- Wilson framework. 
On the other hand, when interactions are strong enough 
compared with kinetic energy, an infinite interaction hmit 
can be a good starting point. The presence of such a large 
energy scale gives rise to a constraint in dynamics of elec- 
trons. In addition, interesting physics now arises from 
the kinetic-energy contribution in the restricted Hilbert 
space, and "non-local" order parameters or more care- 
fully, link variables instead of on-site ones in lattice mod- 
els appear as important low energy collective degrees of 
freedom. These link variables can be formulated as gauge 
fields, and gauge theory arises naturally for dynamics of 
strongly correlated electrons. 

In the present study we revisit a doped Mott insulator 
problem based on one possible gauge theory approach of 
the slave-fermion representation. The slave-fermion rep- 
resentation relies on the fact that a spin degree of freedom 
of an electron is represented by a bosonic matter, thus it 
is advantageous in describing magnetic ordering of local- 
ized moments via condensation of bosonic spinons while 
it has difficulty in connection with Fermi liquid. On the 
other hand, the slave-boson approach resorts to the fact 
that a charge degree of freedom of an electron is repre- 
sented by a bosonic matter, thus the slave-boson gauge 



theory has its clear connection with the Fermi liquid 
theory via condensation of bosonic holons while it docs 
not have a "natural" description for ordering of local- 
ized magnetic moments, compared to the slave-fermion 
approach. [Ij An immediate question is about statistics of 
spinons. At present, there is no clear connection between 
the fermionic and bosonic descriptions for spin degrees of 
freedom. 2] 

The problem to determine which of the representa- 
tions is available seems to be associated with the nature 
of an undoped parent Mott insulating phase. Although 
there is no obvious way to classify undoped Mott insu- 
lating phases, it is natural to consider symmetric and 
symmetry-broken Mott insulators. Q Symmetric Mott 
insulators are usually called spin liquids, where both 
fermionic Q and bosonic Q spin descriptions are avail- 
able. Symmetry-broken Mott insulators can be discrim- 
inated further, depending on their symmetry breaking 
patterns. An example of translational symmetry-broken 
Mott insulators is a Mott insulator with a valence bond 
solid order, and that of rotational symmetry-broken ones 
is an antiferromagnetic Mott insulator. Bond-operator 
formalism!^ may be useful for the description of the 
valence-bond-solid Mott insulator while the antiferro- 
magnetic Mott insulator will be described well by the 
Schwinger-boson representation[^ . 

In this paper we focus on doping to an antiferromag- 
netic Mott insulator based on the slave-fermion approach. 
In early days single-hole dynamics in the antifcrromag- 
netically correlated spin background was one of the main 
interests, [s] Hole doping frustrates an antiferromagnetic 
spin configuration. Interactions with such antiferromag- 
netic fluctuations give rise to a self-energy correction of a 
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single-hole, resulting in the fact that it resides in four di- 
agonal momentum points of (±7r/2, ±7r/2). Mathemati- 
cally speaking, the hopping term in the slave-fermion rep- 
resentation of the t-J model, usually called the Shraiman- 
Siggia term in the continuum approximation, plays a cru- 
cial role for the single-hole dynamics. In the present pa- 
per the role of this term is critically reinvestigated, and 
we find that it plays an important role for dynamics of 
both holes and spin fluctuations. 

Recently, this problem was revisited by several au- 
thors, claiming that the slave-fermion approach allows 
an anomalous metallic phase with short-range anti- 
ferromagnetic correlations, where bosonic spinous are 
gapped away from the antiferromagnetic Mott insulat- 
ing phase. 0, [lO| An effective field theory for the doped 
antiferromagnetic Mott insulator was argued to be a non- 
relativistic fermion (holon) U(l) gauge theory for charge 
dynamics (of gapless charged fermions around the four 
diagonal points) and a relativistic boson (spinon) U(l) 
gauge theory for spin dynamics (the CP^ gauge the- 
ory of the 0(3) nonlinear a model for antiferromag- 
netic spin fluctuations), allowing an anomalous metal 
phase when bosonic spinous are gapped. pH The pres- 
ence of such an anomalous metallic phase is attributed to 
deconfinement[12, 13, 14] of compact U(l) gauge theory, 
where the presence of gapless fermionic matters is ex- 
pected to suppress monopole (skvrmion[l5j) excitations 
in two space and one time dimensions [(2-1- 1)D] although 
the pure compact U(l) gauge theory does not allow such 
a deconflnement phase in (2 -I- 1)I?{1Q]. 

From a microscopic point of view, the emergent U(l) 
gauge structure can be allowed when hopping of spinous 
and holons between next nearest neighbor sites, i.e., hop- 
ping between same sublattices in the square lattice has 
its nonzero vacuum expectation amplitude while hop- 
ping between different sublattices (nearest neighbor sites) 
vanishes in the presence of short-range antiferromagnetic 
correlations (spinon singlet-pairing order). In the present 
study we reexamine the spinon-holon exchange hopping 
term carefully, and reveal that the emergent gauge struc- 
ture is Z2 instead of U(l) owing to the fact that hopping 
of doped holes frustrates the coUinear antiferromagnetic 
spin configuration and results in a complex magnetic 
structure. Such a complex spin texture is reflected in 
nonzero expectation values of nearest-neighbor hopping 
parameters of spinous and holons. In other words, dy- 
namics of doped holes enhances ferromagnetic spin fluc- 
tuations, reducing the U(l) gauge symmetry, arising in 
the case when only coUinear antiferromagnetic correla- 
tions are considered, down to Z2. 

For definiteness, we first consider the antiferromag- 
netic Heisenberg model based on the Schwinger-boson 
representation. In this study we allow both ferromagnetic 
(exchange hopping of spinous between different sublat- 
tices) and antiferromagnetic (singlet pairing of spinous 
between different sublattices) correlations. In addition, 
emergence of a flux is admitted for the ferromagnetic 
exchange-hopping channel. Based on the saddle-point 



analysis for magnetic ordering, we find; (1) If no flux is 
allowed in the exchange-hopping parameter, the vacuum 
expectation value of the hopping parameter vanishes, 
thus our mean-fleld analysis recovers the conventional 
Schwinger-boson mean-field theory for coUinear antifer- 
romagnetic ordering. [tI The resulting effective field the- 
ory is a relativistic U(l) gauge theory, i.e., the CP^ gauge 
theory of the 0(3) nonlinear a model. [17] On the other 
hand, (2) if 7r-flux is allowed in the hopping channel, 
the hopping parameter between the different sublattices 
has its nonzero vacuum expectation value because the 
presence of 7r-flux enhances ferromagnetic correlations. 
As a result, the U(l) gauge symmetry is down to Z2 
since both spin-singlet pairing and exchange hopping ex- 
ist. The presence of both antiferro- and ferro- magnetic 
contributions causes the magnetic ordering structure to 
be complicated. We discuss physics of this possible com- 
plex magnetic structure. 



Next, we study a doped antiferromagnetic Mott in- 
sulator based on the U(l) slave-fermion representation 
of the t-J model. As discussed before, the main point is 
how to take the hopping-t term into account. Performing 
the similar mean-field analysis as the above, we find that 
the hopping parameter does not vanish even in the uni- 
form flux case owing to the fact that mobile doped holes 
increase ferromagnetic correlations. As a result, a com- 
plex magnetic structure such as a spiral type is expected 
to appear, consistent with doping-induced spiral mag- 
netic ordering of the previous slave-fermion studies. [8| 
Enhancement of ferromagnetic correlations induced by 
doped holes causes a serious result that gapless U(l) 
gauge fluctuations do not exist in the long-wave length 
and low-energy limits even when such bosonic spinous 
become gapped to cause the fact that the spiral-type 
complex magnetic structure disappears away from half 
filling. As a result, this spin-gapped anomalous metallic 
state, where only charge fluctuations are gapless exhibit- 
ing so called spin-charge separation, is described by a 
Z2 gauge theory instead of a U(l) gauge theory. In this 
respect such a deconflnement phase is certainly possible 
for (2 -I- 1)D, considering that even the pure Z2 gauge 
theory without matters allows its deconflnement phase 
in (2 -|- 1)1). [31 We propose this gauge-symmetry reduc- 
tion in the low energy limit as one possible mechanism 
for spin-charge separation in the doped antiferromagnetic 
Mott insulator. 
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II. ANTIFERROMAGNETIC MOTT 
INSULATOR 

A. Schwinger-boson mean-field analysis of the 
antiferromagnetic Heisenberg model 

We consider the quantum antiferromagnetic Heisen- 
berg Hamiltonian 



H 
~J 



(1) 



Generally speaking, one can decompose the spin- 
exchange interaction Si 
follows, 

(i) Neel scheme. 



Sj with Schwinger-bosons as 



(ii) Spiral scheme, 



Si ■ Sj — 



S/2 



(2) 



(3) 



where Ai 



biibji)/2 and xlj = (^^V^T + 
bl^bji)/2 are associated with antiferromagnetic and fer- 
romagnetic correlations, respectively. 

In order to clearly show physical meaning of Aij and 
Xij, it is convenient to consider classical spins. If we 
choose two nearest neighbor spins 



/2S 



'2S 



cos(6i/2) 
sin(6i) 



(4) 



where 



are spherical coordinates of spin j, then 

^ Ssin{0/2) (5) 
Xy- =5e'*cos(0/2). (6) 



For a ferromagnetic state {6 = 0), we have A = and 
Ixl = S while in the Neel state of = tt, A = S* and 
X = result. In this respect A and x represent antifer- 
romagnetic and ferromagnetic correlations, respectively. 
Furthermore, for classical antiferromagnetic spins, when 
spins are non-coUinear (i.e. 9 n) where neighboring 
spins have their finite overlap, x becomes nonzero so 
that X describes canting and spiral. One may expect 
that competition between antiferromagnetic A and fer- 
romagnetic X correlations will give rise to non-coUinear 
ordering in the ground state of antiferromagnctically cor- 
related spins. 

Inserting Eqs. (5) and (6) into Eqs. (2) and (3), we 
get 



{Si ■ Sj)Neel = S^ COS 9 
(Si ■ Sj)spiral ~ S^ COS 9 — S/2. 



(7) 

(8) 



ground state energy. In the following we focus on the 
Spiral scheme. 

Taking into account the constraint of ^lo-^io- — 2S 
at each site via a Lagrange multiplier field A^, an effec- 
tive Hamiltonian in the Schwinger-boson representation 
of the antiferromagnetic Heisenberg model can be written 
as follows 



H 



J = -^{A*jCrh„bja + h.c 
+ ^A,(6L6..-25) + ^(|A 



IXijl 



h.c 



(9) 



It is important to notice that the sign of the spinon- 
hopping term is positive, implying that spinon-hopping 
contributions are energetically unfavorable, thus its am- 
plitude will vanish in the uniform-flux case, as will be 
shown below. Physically speaking, ferromagnetic cor- 
relations frustrating the antiferromagnctically correlated 
background are prohibited in the naive mean-field anal- 
ysis. However, introduction of 7r-flux weakens such frus- 
tration effects, ferromagnetic correlations contributing to 
physics of antiferromagnetism. 

For a mean-field analysis we replace Ay and Xij with 
their saddle-point values. In the flux ansatz we parame- 
terize Ay and Xij with 



A,- 



X^j = xe 



(10) 



where 9ij = ±9 if ij along or against an arrow in Fig. 
1. Note that a flux is allowed only in the ferromagnetic 
hopping channel. Based on this mean-field ansatz with 
Xi = A, we obtain the following effective Hamiltonian in 
the momentum space. 



"7 = --j-^^^^'^^'^^^i^'f 



h.c.) + Xj^bl^bka 



'Zx 



Ik COS 9 
—i(f>k sint 



ka 

i(j)k sin 6* 

-7A; COS ( 



bka- 

bk+Qc 



-2NSX 



-TVA^ 



'-Nx^ 



(11) 



>k e 



10 



IJ 



It is obvious that the Spiral scheme is more favorite than 
the Neel scheme in the classical limit since it has lower 



FIG. 1: Flux ansatz 
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with 7fe = ^/zY.jeM'''' + e''"") and (j)k = 
1/2; ^^g^(e*'^"= — e*'^«). z is a lattice-coordination num- 
ber, here z = 4 for the square lattice. The ' symbol 
in the momentum summation represents to perform the 
summation in the folded Brillouin zone. 

Performing the Bogoliubov transformation, we can di- 
agonalize the effective Hamiltonian as follows, 

fe,s=± 

- 2N\{S+]^)+'-N/:^^-'-Nx\ (12) 

where the quasiparticle spectrum is 

= ^{\ + szxQk?-{zMk/2y (13) 

with 0fe = (7fe cos 6Y + {(pk sin 9Y in the flux ansatz. 
Then, the resulting mean-field free energy is given by 

ks 

-2NJ\{S+\)+'-NJI\''-'-NJx^. (14) 

Minimizing Fmf with respect to A, x, and A, we ob- 
tain self-consistent mean-field equations 

5+1/2. l^A±i2^(„, + 1/2) (15) 

a4i:^k+v2) (16) 

ks ^ 

1 ^ 2zs9fc(A + sxQk) , .^„. 
^=^2^ -s (^fe + l/2) (17) 

with A — > zA/2 and x ~^ ^x- 

Flux states of our interest are uniform-flux [d = 0) and 
TT-flux {Q = 7r/4) phases since they do not break time- 
reversal symmetry. For the uniform phase, our numeri- 
cal analysis shows that the lowest free energy is achieved 
by A = 1.1580 and x = 0- One can easily check that 
X = is a solution for the above self-consistent equations 
in the uniform-flux case since the right-hand-side of Eq. 
(17) is odd for momentum fc if x = 0, thus it vanishes 
after the momentum integration, consistent with the left- 
hand-side of Eq. (17). In addition, this result is consis- 
tent with our physical picture associated with the sign of 
the ferromagnetic hopping term. Schwinger bosons are 
condensed to the k* = (0, 0) state, which corresponds to 
the collinear antiferromagnetic order. Staggered magne- 
tization is found to be mg = 0.3034, which is the same 
as the Auerbach's result [7]. The ground state energy per 
bond is given by 

Eo = -A^ = -0.3352. (18) 



+ 2 - 2 




-3 -a -i 1 2 3 -3-2-1012 



FIG. 2: (Color online) oj^'^ and tJj. ^ in the vr-flux phase when 
A = A - y2x/2 with A = 1.2021 and x = 0.2349, where there 
exists a contour for the quasiparticle spectrum to cross from 
positiveness to negativeness in ijJ^^. 

For the n-flux phase, there are two branches of bosons, 
tj^l^-boson and w^-boson. Spinon condensation arises in 
each branch when 

(1) A + V2x/2. At the wave vector k* = (0, 0) we 
obtain w^, > but uj^, — 0. The self-consistent equation 
(fTT]) deduces x = corresponding to the lowest energy, 
which is the same as the uniform-flux case, and 7r-flux 
doesn't give any particular effect. Neel ordering occurs 
only for the w^T branch. 

(2) \ — A ~ \/2x/2. CLij!"-branch bosons are condensed 
to the k* = (0, 0) state, which implies that collinear 
antiferromagnetic ordering is also realized within uj'^ 
bosons. On the other hand, if we examine the dis- 
persion of low lying excitations of uj^ bosons, we find 

uj^ w V2A/2 (|fcp - 4V2x/A)^^^ around the zero point. 
That is, there is a ring divergence rather than points pro- 
vided that X 7^ (see Fig. Thus, for the branch 
the bosons with < 4\/2x/A are all going into the 
condensate portion of the boson state. In other words, 
the density of uj^ bosons stays at |A:p < 4\/2x/^ '^^ mo- 
mentum space, thus many kinds of non-coUinear mag- 
netic ordering emerge. Actually, we obtain a non-zero 
X = 0.2349 and A = 1.2021, solving Eq.!!?]) numerically. 
Moreover, comparing with the ground state energy of the 
uniform-flux phase Eq. , this complex magnetic state 
is found to have a lower energy indeed, 

= -A^ + x^ = -0.3578. (19) 

We can extract the collinear staggered magnetization 
from the branch, and find tuq = 0.1043. This value 
is much smaller than that of the uniform-fiux case owing 
to enhancement of ferromagnetic correlations. 

B. Discussion 

In the previous section we have seen that ferromag- 
netic correlations do not contribute to physics of collinear 
antiferromagnetism in the conventional Schwinger-boson 
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mean-field approximation, i.e., uniform-flux ansatz. On 
the other hand, the tt-Aux phase turns out to be more 
stable than the uniform-flux phase, and the emergence 
of TT-flux enhances ferromagnetic correlations, causing a 
complex magnetic state, where there exist various order- 
ing wave vectors. 



III. DOPED ANTIFERROMAGNETIC MOTT 
INSULATOR 

A. A slave-fermion mean-field theory of the t-J 
model 

1. A slave-fermion effective action 



Nature of the magnetically ordered state in the tt-Aux 
phase is an important question. The present semiclassi- 
cal approach does not seem to be sufficient for determin- 
ing the resulting ground state. This situation reminds us 
of "order-by-disorder" ,(T^ where quantum fluctuations 
split degeneracy and determine its ground state. One 
possible scenario in the square lattice is that magnetic 
ordering is still to show collinear antiferromagnetism de- 
spite the presence of such "degenerate" states, i.e., many 
A;*s corresponding to ordering vectors, since magnetism 
measured from the degenerate states will be cancelled out 
after momentum integration if the condensation ampli- 
tude at each momentum point is assumed to be the same 
as each other. These contributions reduce the strength of 
the collinear antiferromagnetism, as seen in the previous 
calculation. This is an important effect of ferromagnetic 
fluctuations in the antiferromagnetic ground state. 



Existence of the complex spin texture prohibits emer- 
gence of gapless U(l) gauge fluctuations owing to nonzero 
vacuum expectation values of both spinon-hopping and 
spinon-pairing order parameters. Instead, Z2 gauge the- 
ory appears for possible non-magnetic phases resulting 
from this complicated spin pattern. Q Then, deconflned 
spinon excitations are expected to appear at high ener- 
gies beyond their excitation gap. Actually, dynamic spin- 
excitation spectra have shown that there is large spectral 
weight at high energies be yon d the spin- wave approxima- 
tion in undoped cuprates. 20] In the present context this 
unidentified spectral weight may be identified with de- 
confined spinon continuum allowed by Z2 gauge theory. 



However, it should be noted that the Z2 gauge struc- 
ture is allowed only in the low energy limit, more pre- 
cisely, T < Ty^, where the ferromagnetic exchange- 
hopping parameter has its nonzero expectation value be- 
low T^. Then, U(l) gauge fluctuations appear above this 
temperature. The Z2-U(l) crossover makes deconfine- 
ment of spinous not clear. This suggests to measure dy- 
namic spin spectra carefully at high energies as a func- 
tion of temperature although it is beyond the scope of 
the present mean-field study to predict this temperature 
accurately. 



We start from the t-J model for doped Mott insulators 

ij ij 

(20) 

Based on the U(l) slave-fermion representation 

Cia = l/jlbtcr (21) 

with the single occupancy constraint ^L^jo- + V"! V'i = 
25 where S* is a value of spin, here S = 1/2, one can 
rewrite the spin-exchange term and electron hopping 
term in terms of fermionic holon and bosonic spinon op- 
erators as follows 

J '^{S^ ■ Sj - ^n,nj) 
ij 

ij ij 

-tJ2(^lc,. + H.c.) ^ tY,{xU% + H.c.) 

ij ij 
ij ij 

(22) 

where the bond operators 

a(3 (J 

(23) 

describe antiferromagnetic, ferromagnetic correlations, 
and hopping of doped holes, respectively. Comparing 
with the Heisenberg Hamiltonian [Eq. ([5])], one can see 
— ^Xij and J l^ij correspond to Jxij and JAy , respec- 
tively. 

Using the above expressions, the t-J model in the elec- 
tron language is mapped onto exactly the same model but 
with fractionalized particles such as spinons and holons 
in the following way 

^e// = E ^U^r - Ai)&.. -tY^{hlx%h„ + H.c) 

i ij 

ij ij 

i ij ij 

E ^Mab^a + " 2S), (24) 



where is a Lagrange multiplier field to impose the sin- 
gle occupancy constraint. 

It is important to notice that ferromagnetic exchange 
correlations in the spin-exchange term is not allowed in 
the hole-doped case. This implies that at half-filling, 
this effective Lagrangian reduces to the conventional 
Schwinger-boson theory instead of the "flux" Schwinger- 
boson theory since the absence of holons results in 
Xij — 0, causing ferromagnetic correlations to vanish, 



Xi 



0. The reason why we do not introduce the fer- 



romagnetic correlation term from the Heisenberg term is 
that when holes are doped, the presence of charge contri- 
bution, — jniUj, does not allow the simple ground-state- 
energy analysis even for the case of classical spins, per- 
formed in the study of the Heisenberg model, thus it is 
not obvious that introduction of such contributions low- 
ers the ground state energy. Here, we focus on the hop- 
ping term of the t-J model instead, giving rise to ferro- 
magnetic spin correlations induced by doped holes. This 
doping-induced ferromagnetic-correlation term will play 
an important role in the structure of magnetic ordering. 
As seen from its sign, opposite to the half-filled case, it 
will contribute even in the uniform-flux case, and modify 
the collinear antiferromagnetism drastically. 

In the following we will consider the 7r-flux phase al- 
though we comment on the uniform-flux case briefly. 
There are two reasons for considering such a flux phase. 
The TT-flux phase is time-reversal symmetry-preserving, 
thus physically allowed. More importantly, the tt-Aux 
ansatz gives rise to four hole pockets, consistent with 
the previous studies Q, where the spinon-holon exchange 
hopping term (Shraiman-Siggia term) is taken into ac- 
count in the context of the perturbation theory, reveal- 
ing that holons reside in the four diagonal points. Al- 
though we do not know the connection between our flux- 
phase result and the previous perturbation theories, the 
TT-flux ansatz reproduces such hole pockets. [2l[ On the 
other hand, the uniform-flux phase does not allow the 
hole pockets. 




FIG. 3: (Color online) Energy dispersion of fermionic holes, 
where Min[^f] = results at = (±7r/2, ±7r/2). 



spin order in the presence of hole doping. Together with 
the TT-flux ansatz 



ip i) tp fp _ A -i 

□ □ 



(26) 



the free energy can be decomposed into the charge 
(fermion) and spin (boson) sectors, 

(1) Fermion sector 



2. n-flux ansatz 

Four order parameters are assumed to be as follows for 
the mean-fleld analysis 



Xij 



y Af,, iXi 



-A, 
Xij 



(25) 



decoupling the slave-fermion effective Lagrangian, where 
the spin-singlet pairing order parameter and Lagrange 
multiplier field are assumed to be uniform and real while 
the hopping parameters of spinous and holons have their 
phase contributions. A nonzero value of Af, corresponds 
to short-range antiferromagnetic correlation as discussed 
before, and a nonzero Xb iXi') implies hole (spinon) mo- 
bility or ferromagnetic correlation, which deforms the 



^'A-~^EK2cosh(^)} 

ks 

/ 

+2txi^ V(7fe cos 6)2 + (v3fcsine)2 

= - ^ E [^^{ 2 cosh ( I [2tx^ V(7fcC0se)2 + ((^fc sin 9)2 

k 

+A])} +ln{2cosh(^[-2ixv.V(7fcCOse)2 + (^fcsine)2 
I }] + 2txri, VhkCOsSy + iipk sin 9)2 (27) 



(2) Boson sector 



7 



ksl 

= I 51 sinh(^{ (2tx6 V(7fc cos 9)2 + (^^ sin 6)2 

A; 

-M + A)'-(4JA,7.)'}'^')} 
+ ln{2 sinh(| { (-2txb v/(7fc cos 6)2 + {ip^ sin 9)2 
-M + A)^-(4JA,7.)'}'^')} 



(28) 



where 



7fc = cos kx + cos fcy, </3fe = cos kx — cos fcj,, 46 = tt. 

(29) 

Noted from the energy dispersion = 
{"fk COS 6)2 + ((^j, sin 6)2 of fermions (Fig. [3]), 
one can see holes are pocketed around the diagonal 
points, (±7r/2,±7r/2). 



3. Self-consistent equations 

In this study we focus on the properties at zero tem- 
perature. Minimizing the free energy with respect to A;,, 
Xbj Xi>^ A and /i, we obtain the following self-consistent 
equations, 

(1) Fermion sector 



5 = 



1 



2Nl 



E 1 



1 



Xb 



(30) 
(31) 



Above equations clearly show that holons are most 
densely located within the Dirac cones at fc = 

(±7r/2, ±7r/2). For small doping we have S — ^ (si^) 

and Xb = —^^S^^^, as derived in appendix. 
(2) Boson sector 



saddle-point equations are reduced to the conventional 
Schwinger-boson mean-field equations in the zero-doping 
limit {S — > 0) since Xb vanishes from Eq. (31) and XV' 
disappears accordingly from the summation for s — ± 
in Eq. (32). Hole doping effects reduce the amplitude 
of an effective spin from 25" to 25 — (5, as shown in Eq. 
(34), and increase ferromagnetic correlations, as shown 
in Eq. (32). Accordingly, antiferromagnetic correlations 
become weaken, seen from Eq. (33). 

In the mean-field treatment magnetic ordering is de- 
scribed by condensation of Schwinger bosons. Here, we 

introduce XbC' ^bc ^^'^ ^bc to describe the contribu- 
tion from the E\{k*) — states, which are taken into 
account separately: 



- E 



-sC{\ + '2stxbe,) 



mk) 

'Xbc 2^ = XV- (35) 

k* ,s — ± 



= {2S+l)-5. (37) 

If E'^ik*) = and Ek',s=±riBc{2stxb^t + A)/A is fi- 
nite, spinous are condensed to the E'^{k*) = states, and 
the system possesses a magnetic long-range order. The 
ordering wave vector k* is determined by E^{k*) = 0, 
and especially, k* = (0, 0) corresponds to long-range 
anti-ferromagnetic ordering. For definiteness, we have 
checked that our treatment of these mean-field equa- 
tions recovers the uniform-flux Schwinger-boson mean- 
field theory correctly. At half-filling the spinon spectrum 
of El{k) = Etik) is gapless at fc = (0,0) which en- 
sures their Bose condensation. Solving Eq. (jST]) . we see 
that the spontaneous staggered magnetization is given 
by {S^) — riBC ~ 0.3, recovering the half-filled result 
completely. 



1 ^Zftflil^,,, ,32, 



2A^L 



k,s=± 



Em 



-E 



^J^bll 



-Ah 



k,s=± 



(33) 



1 ' [istxhit + A 

sri:^^|r^ = (2s+i)-* (34) 

with the spinon-quasiparticle spectrum E'^{k) — 
Y (2stx6?fc + A)^ — (4JA67fc)^. It is clear that these 



4- Phase diagram 

Now, let's consider how the antiferromagnetic long- 
range order evolves with doping. Due to the 7r-flux order, 
there are two types of bosons, with energy —E'^{k) 
and B- with energy —E^_{k). Away from half- filling 
holes are constrained within the Dirac cones, thus Eq. 
([3T|l gives Xb < 0, which deduces B\_ has lower energy 
than Bt at finite doping, i.e., -E\{k) < -Et{k). When 
A 8 JAfc + 2^/2txb is satisfied, we have E^_ (0) for 
B^ bosons, thus their Bose condensation results at the 
wave vector k* = (0,0), implying an antiferromagnetic 
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aM 



J/t=3 



0.2 0.3 

6 



FIG. 4: (Color online) Doping dependence of A;, and ubc 
with t/J = 3.j22] In order to clearly show the tendency, 
and Ubc are scaled with the values at half-filling. The spin- 
gapped anomalous metallic phase is marked as aM. 



order. For bosons within the same value of A, one 
can easily check that 1/E'^{k*) shows divergence along 
a ring k* G {k\E']^{k) ~ 0}. In other words, the den- 
sity of -B+ spinons resides at large number of k* states, 
thus many kinds of incommensurate magnetic ordering 
emerge and such a complex spin configuration deviates 
from the coUinear antiferromagnetic ordering, similar to 
the half-filling case in the mean-field approximation. 

With increasing doping and t, electrons become more 
active, implying that ferromagnetic correlations de- 
scribed by Xb and Xip should also increase. However, 
the antiferromagnetic order represented by the conden- 
sation amplitude ngc becomes more suppressed with 
large i/J, which is confirmed by the numerical results. 
Furthermore, by numerically solving the self-consistent 
equations for given value of t/J = 3, Fig [4] exhibits that 
antiferromagnetic ordering disappears completely above 
Sc ~ 0.25 22] while the spin-singlet order persists up 
to much larger hole concentration (Af, ^ 0) away from 
the antiferromagnetic order. We identify this intermedi- 
ate non-magnetic spin-gapped phase with an anomalous 
metal, as will be discussed in the next section. 

Compared to the uniform-flux case although not shown 
in the present paper, we have qualitatively the same 
phase diagram. This originates from the fact that the 
sign of the ferromagnetic correlation term is negative 
opposite to the Heisenberg model, thus such contribu- 
tions are energetically favorable even in the uniform-flux 
phase. However, the presence of the 7r-flux order en- 
hances ferromagnetic correlations as seen in the study of 
the Heisenberg model. As a result, the coUinear antifer- 
romagnetic order is killed more rapidly via hole doping 
in the vr-flux phase. We have also examined the role 
of frustrated hopping effects, i.e., next-nearest-neighbor 
hopping terms. Such frustrated hopping effects modify 
bare dispersions for spinons and holons. However, we do 
not find any qualitative changes in the phase diagram 
since the presence of 7r-flux dominates over the band- 
modification effect. 



B. Spin-gapped anomalous metal 

As shown in the phase diagram of Fig. 4, hole doping 
kills the antiferromagnetic order and results in a para- 
magnetic state, where bosonic spinons are gapped but 
short-range antiferromagnetic correlations still remain. 
We identify this phase with an anomalous spin-gapped 
metal since charge excitations carried by holons are gap- 
less, but spin fluctuations are gapped, thus exhibiting 
spin-charge separation. 

To understand physics of this anomalous metallic 
phase, it is necessary to derive an effective field theory. 
Considering holons first, one can expand the holon field 
as 



+ e 



ik2X 



-ik2X 



■02(2;) 



(38) 



around the diagonal points fci — (7r/2,7r/2) and ^2 = 
(7r/2, — 7r/2), where the holon spectrum is linearized as 



E]^ = 2tx^\ cos^ kx + cos2 ky sa v^J kl + kl (39) 



with its velocity 

by 



"itxip- The holon spinor is given 



V'i(a;) 



1p2o{x) 
i>2eix) 



,(40) 



where e and o represent even and odd sites, respectively. 
Such two-component spinors can be combined to form a 
single four-component spinor as 



Tpix) 



■>pi{x) 

1p2ix) 



I i^le(x) 
•iploix) 
ll^2o{x) 

\ 1p2eix) 



(41) 



An important point is that U(l) gauge fiuctuations be- 
come gapped due to condensation of a "charge" 2 scalar 
field, here the ferromagnetic hopping parameter, reduced 
to gapped Z2 gauge fields, as discussed previously. In 
other words, hole doping causes magnetic frustration ef- 
fects, giving rise to ferromagnetic correlations. This re- 
duces the emergent U(l) gauge structure down to Z2. As 
a result, we find the following effective field theory for 
the anomalous Z2 spin-gapped metal phase 



(42) 



where Dirac gamma matrices satisfying the Clifford al- 
gebra [7,1,7,.]+ = 26^^ are given 



70 



( era 






f 0-2 






I 




, 71 = 


V 


-0-2 1 


, 72 = ( 



CTl 







-(Tl 



and fir is a renormalized chemical potential associated 
with hole doping. is an electromagnetic field. Note 
that gapped spinon excitations are ignored for low energy 
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physics of this metaUic phase. Then, Eq. (42) is our 
effective field theory for the anomalous metallic state. 

An immediate question is about the transport prop- 
erty in this metallic phase. Since there is no scatter- 
ing mechanism with gapless gauge fluctuations in the 
above effective field theory, Fermi liquid behavior is 
expected to appear in the low temperature limit, i.e., 
5p{T) = p{T) - p{T = 0) ~ r2 with resistivity p{T). 
Although this is true in the present description actually, 
such a Fermi liquid behavior in transport is not a full 
story. 

As discussed in the Heisenberg model, there is an en- 
ergy scale associated with the crossover behavior from 
Z2 to U(l) since ferromagnetic contributions are nonzero 
below this temperature. Although we cannot determine 
this temperature satisfactorily in our mean-field analy- 
sis, gapless U(f) gauge fluctuations should arise above 
this crossover temperature. Then, scattering with such 
gapless gauge excitations would lead to non-Fermi liquid 
physics in this regime. Frankly speaking, there are var- 
ious regimes at finite temperatures, associated with not 
only the Z2-U(l) crossover temperature but also so called 
the spin-gap temperature T^, where bosonic spinous are 
effectively critical above this temperature. Tf, should not 
be confused with Ta, where spinon-singlet pairing order 
appears below this temperature, much larger than the 
spin-gap temperature Ti,. In this respect we should con- 
sider both energy scales, and classify various finite tem- 
perature regimes. In the present paper we discuss only 
low energy physics of this metallic phase. 



IV. DISCUSSION AND SUMMARY 

A. Coherent electron quasiparticle excitations near 
the diagonal points 



In the low temperature regime, more precisely in the 
zero temperature limit where we have studied in this pa- 
per, the situation becomes worse because even gapless 
U(I) gauge fluctuations do not exist owing to the gauge- 
symmetry reduction. However, we argue that such an 
electron excitation may still appear as a resonance state 
between a spinon and holon. To see this resonance state, 
we revisit the hopping-t term in the slave-fermion repre- 
sentation of the t-J model. We rewrite the hopping term 
as follows 

-tY^^^lhlha + H.c.) ^ tY,{clc,, + H.c.) 

ij ij 
ij ij 
ij ij 

+tY,{^,]x%^, + H.c.). (43) 

ij 

It should be noted that the electron operator Cia is 
different from the original one since it is a Hubbard- 
Stratonovich field although it is a fermion. Rather, it 
should be understood as a collective fermionic excitation 
in the spinon-holon liquid. We identify this fermion field 
as an electron, which is expected to be seen in the low 
energy regime. Note that hopping of such low energy 
electrons is frustrated, viewed from its sign. 

The electron-spinon-holon exchange term reminds us 
of the slave-fermion representation of the Anderson lat- 
tice model, where localized spins are represented in slave- 
fermions.i27] To find a self-energy correction from this 
"hybridization" coupling seems to be an interesting and 
important problem since the spin-charge separation sce- 
nario may explain the coherent nodal electron excita- 
tions. 



The present slave-fermion mean-field description can- 
not explain the emergence of sharply defined electron 
quasiparticles near the diagonal points and their Fermi 
liquid behavior because electrons decay into spinous and 
holons in the Z2 anomalous metal phase. [24] This is the 
main criticism in this kind of decomposition approach if 
one tries to connect the spin-charge separation scenario 
with pseudogap physics of high Tc cuprates. However, 
if we consider the regime of T > T-^, the situation is 
not perfectly clear owing to the presence of gapless gauge 
fluctuations beyond the mean-field description since such 
gauge fluctuations cause attractive interactions between 
holons and spinous, binding them, which may allow a 
resonance state corresponding to an electron-like excita- 
tion although it differs from confinement. Actually, this 
was a little bit demonstrated, where noncompact gauge 
fiuctuations give rise to attractive interactions between 
"particle" and "hole" (holon and spinon in the electron 
spectrum), increasing coherence. [26] We should confess 
that this argument is far from consensus even in the qual- 
itative point of view. 



B. d-wave superconductivity 

Another interesting subject is how to describe d-wave 
superconductivity in this context. Considering the slave- 
boson case, one may regard this study as a straight- 
forward job. In the slave-boson theory condensation of 
bosonic holons results in superconductivity in the pres- 
ence of spinon singlet pairing since they form Cooper 
pairs. However, this structure cannot be applied to the 
slave-fermion context. Fist of all, charged holons are 
fermionic, thus their condensation is forbidden. Then, 
the next choice will be pairing of such fermionic holons. 
Actually, such pairing interactions are usually attributed 
to the presence of gapless U(I) gauge fiuctuations since 
they cause attractive interactions between holons living 
in different sublattices. What we would like to do is 
to construct a mean-field theory for such pairing order 
from the microscopic model itself. However, if we start 
from the t-J model, such an attractive pairing interac- 
tion channel is difficult to find since the exchange term 
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is represented only in the bosonic-spinon language ow- 
ing to the single occupancy constraint, analogous to the 
fact that the J term is also expressed in terms of only 
fermionic spinous in the slave-boson context. 

To avoid this difficulty, we start from so called BCS- 
Hubbard modelfli] 

H = -tY,{clcj, + H.c.) - ^(Ate„^c,„c,^ + H.c.) 

ij ij 

+ i^|Ayf + i7^n,tn»i (44) 

where the Hubbard-Stratonovich transformation for the 
Cooper channel is performed before the constraint is ap- 
plied. This model has its own interest in the fact that we 
can investigate what happens in BCS superconductivity, 
increasing local interactions U . Our strategy is that such 
strong correlation effects are introduced using the slave- 
fermion representation in the large-?/ limit. Then, the 
pairing term is expressed in the slave-fermion represen- 
tation as 

ij 

^ 5^(A|, A^A^ + H.c.) - J^iAl^^^e^M, + H.c.) 

ij ij 

-E(44^1^]+^-^-), (45) 

ij 

where pairing interactions are renormalized due to Hub- 
bard interactions. Replacing A|^A^j and Al^Afj with 

A^^ and Afj, respectively, we find the slave-fermion ef- 
fective Lagrangian for the BCS-Hubbard model in the 
large- [/ limit 

Leff = E blidr - /i)6.. -tJ2{blx1,b,. + H.c.) 

i ij 

- '^{A'^^^eapbiabjp + H.C.) 

ij 

i ij 

-Y,{A%4^]+H.c.) 

ij 

+ Y.^^jr^+H.c. + l\A.,\')+tJ2ixtx'^,+H.c.) 



the absence of doped holes, causing A* = while A 
disappears at large doping, killing A^^ . 



C. Summary 



In this paper we have considered dynamics of doped 
holes in the antiferromagnetically correlated spin back- 
ground. The point is that doped holes frustrate the 
collinear antiferromagnetic spin configuration, resulting 
in more complex spiral-like spin patterns. This reduces 
the U(l) gauge symmetry down to Z2, thus the effective 
field theory appears to be a Z2 gauge theory instead of 
U(l) in the low energy limit. As a result, deconfinement 
of spinous and holons is naturally allowed in (2 + l)D. 
When spinous become gapped via further hole doping, an 
anomalous metallic phase results with spin-charge sepa- 
ration. We argued that there is a crossover energy scale 
Ty. from Z2 to U(l) above this temperature, where trans- 
port should be more carefully studied in the presence of 
such gauge interactions. We also discussed how coher- 
ent electron excitations near the diagonal points can ap- 
pear in the slave-fermion context. Lastly, we have spec- 
ulated superconductivity within the slave-fermion repre- 
sentation. 



APPENDIX: ANALYTIC EXPRESSIONS FOR 
FERMION SELF-CONSISTENT EQUATIONS 
NEAR ZERO DOPING 



Although investigation of this effective Lagrangian is 
beyond the scope of this paper, one can read an impor- 
tant aspect of this Lagrangian immediately. An observed 
dome-shaped superconductivity line is expected to ap- 
pear since A^- vanishes in the zero doping limit owing 



In this appendix we show the analytic expression of Xb 
in the low doping limit. Minimizing the free energy Eq. 
(27) for x^, we find an equation for Xh and obtain the 
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following expression 



Xb 



E 



Cfc<-(P+>)/2tX^ 



+ 



E a 



«fc<(M+A)/2txv, 
1 



E i 

?fc>(A'+A)/2tx^ C*'>-(M+A)/2tx^ 

-( E 4^+ E 4^' 

(A<+A)/2tx^<C;f'<-(M+A)/2tx^ ?^>-(M+A)/2tx^ 

■ E E 

?fc<(A'+A)/2txv. (M+A)/2txv.<e;f<-(M+A)/2tx^ 



E ^: 



C<(A'+A)/2tx^ 



1 



1 



E 4^- E 4^ 

5;f>-(M+A)/2txv, C^<(M+A)/2tx^ 



E 



5*'>-(M+A)/2txv 



Y [E E 



4iV. 



2Nl V 2 if 



2''^^/^^<-(M+A)/2tx. ^SttP^'^ 



2A^L ^ 



1 

2^ 



1 /A* + A\3 



2A^L 



with = ^cos^ /ca; + cos^ ky. 



(A.l) 



In the same way we find the following expression from 
the variation of the free energy for the chemical potential 



-5 = 



1 r 



E 



E 



C^'>-(M+A)/2txv, C;|'<-(p+A)/2tx^ 

E - E ] 

?;f<(M+A)/2txv, C*>(M+A)/2txv, 



( 



«fc >-(A'+A)/2tx^ (A'+A)/2tx*<«^<-(A'+A)/2tx^ 



- E )+ E 

«^>-(M+A)/2txv. ?fc<(A'+A)/2txv, 

■( E + E 

(M+A)/2txv,<«r<-(A'+A)/2iXY. «t<(A'+A)/2tx^ 



2A^r, 



(p+A)/txv,<C<-(A'+A)/2tx^ 



E 



0<«^<-(M+A)/2tx^ 



1 



— X 4 X 



C^<-(M+A)/2tx^ 



1 /U + 2 
^ X TT . 



As a result, we obtain 



(A.2) 



(A.3) 



In the uniform-flux case we find Xb — (5^ smaller than 
the above in the zero doping limit 
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